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00 ■ Abstract 

Ph, We demonstrate that a certain class of low scale supersymmetric "Nelson- 

D . Barr" type models can solve the strong and supersymmetric CP problems 



while at the same time generating sufficient weak CP violation in the — 
system. In order to prevent one-loop corrections to which violate bounds 
coming from the neutron electric dipole moment (EDM), one needs a scheme 
for the soft supersymmetry breaking parameters which can naturally give suffi- 
cient squark degeneracies and proportionality of trilinear soft supersymmetry- 
breaking parameters to Yukawa couplings. We show that a gauge- mediated 
supersymmetry breaking sector can provide the needed degeneracy and pro- 
portionality, though that proves to be a problem for generic Nelson-Barr mod- 
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els. The workable model we consider here has the Nelson-Barr mass texture 
enforced by a gauge symmetry; one also expects a new U(l) gauge superfield 
with mass in the TeV range. The resulting model is predictive. We pre- 
dict a measureable neutron EDM and the existence of extra vector-like quark 
superfields which can be discovered at the CERN Large Hadron Collider. Be- 
cause the 3x3 Cabbibo-Kobayashi-Maskawa matrix is approximately real, 
the model also predicts a flat unitarity triangle and the absence of substantial 
CP violation in the B system at future B factories. We discuss the general 
issues pertaining to the construction of such a workable model and how they 
lead to the successful strategy. A detailed renormalization group study is then 
used to establish the feasibility of the model considered. 
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I. INTRODUCTION 



The strong CP problem is without question one of the most important problems faced by 
the Standard Model (SM). Its origin lies in the necessity of adding the so-called 9 term to the 
effective QCD Lagrangian due to the contribution of instantons present in the topologically 
nontrivial QCD vacuum 

^eff = -^F^^F^^^ , (1.1) 

where the dual field strength is given by F^^ = \^^uai3F°'^ . Through the anomaly in the axial 
U(l) current of QCD, chiral U(l) transformations lead to shifts in 9, leaving the physical 
combination 9 = 9 — argdetMg, where Mq is the quark mass matrix. Since Ces clearly 
violates CP, it gives a strong interaction contribution to the neutron electric dipole moment 
and leads to the experimental contraint 

9 < 10-^ . (1.2) 

The real problem therefore is one of naturalness or fine-tuning: Why is 9 so incredibly small? 

There are currently three notable classes of possible solutions to this problem: (1) van- 
ishing up quark mass, (2) the axion [^,^ or (3) CP conservation and subsequent spontaneous 
breaking. The first and simplest possibility appears to be disfavored by current algebra rela- 
tions between pseudoscalar meson masses but is still controversial (see e.g. Ref. [^). Of 
these, the most popular is the invisible axion alternative [Q. Here one introduces a global 
chiral U(1)pq (Peccei-Quinn 0) symmetry which is spontaneously broken at a high energy 
scale / and explicitly broken by instantons. The 9 parameter is replaced by a dynamical 
field — a pseudo Goldstone boson of the U(1)pq — whose potential dynamically relaxes 9 
to zero. The advantage of this scheme is that it is simple, generic and has observable con- 
sequences both in terrestrial experiments and in astrophysics and cosmology. Astrophysical 
constraints from axion-induced cooling during stellar evolution ^ and effects on the neutrino 
signal from supernova 1987A give a lower bound, / > 10^° GeV, while a cosmological 
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upper bound of ~ 10^^ GeV is given by the contribution to the universal energy density of 
the vacuum energy associated with U(1)pq breaking as the axion vacuum expectation value 
relaxes to zero [|iy]. On the aesthetic side, one may complain that we are merely replacing 
the 9 fine-tuning problem with another: the smallness of the ratio of the weak scale to 
the U(l)pQ breaking scale ~ 10"'-^"^'^''. Another possible problem is the dependence of the 
solution on a global symmetry, generally not preserved by gravity, so not likely to appear 
from a more fundamental theory. This appears to be a significant problem at least in 



Einstein gravity. However, it has been argued that gravitational violations of global symme- 



tries may be suppressed in certain extensions, including string theories [|12|, where at least 
the universal dilaton-axion is always present. The axion alternative also does not provide 
an explanation of weak CP violation. Here one must assume the Kobayashi-Maskawa (KM) 
origin of CP violating phases. Of course the ultimate test is to detect actually an axion ||13|| . 

In this paper we will focus on the third alternative. That is, we will assume that the fun- 
damental theory of nature preserves CP and that at sub-Planck energies it is spontaneously 
broken. Indeed there is evidence that CP has its origin as a gauge symmetry remnant of 



superstring theories |]T^. In this manner the smallness of 9 reflects the existence of an 



underlying symmetry. Such models were first constructed in the context of Grand Unified 
theories (GUTs) by Nelson and refined by Barr |jl5| and incorporate extra heavy quarks 



which mix with the observed quarks. Then, relying on specific symmetries, one can obtain 
a texture of the full quark mass matrices which guarantee the tree-level vanishing of 9 after 
the spontaneous CP violation (SCPV). After integrating out the heavy fields, the low energy 
quark mass matrices contain the usual KM phase. The generic difficulty with these models 
comes from the need to ensure that large contributions to 9 do not arise at higher loops, 
while at the same time having sufficient weak CP violation from the KM phase. Thus while 
the SCPV approach is conceptually rather simple, it is not so generic and requires careful 
model building. 

Given some of the tantalizing hints of low energy supersymmetry and the plausible 
gauge origin of CP symmetry, it is worthwhile to attempt to construct a supersymmetric 



(SUSY) model with a Nelson-Barr type mechanism for solving the strong CP problem. The 
SCPV feature then also resolves the so-called SUSY phases problem. The latter problem 
was originally described in the context of a minimal supergravity origin of the soft SUSY 



breaking terms [|T^], and is usually worse in a general SUSY breaking scenario. There are 
two extra phases in the universal soft mass parameters, beyond 6km and ^, which give CP 
violating effects in the low energy effective theory. These can be written as effective phases 
in the coefficients A and B of the trilinear and bilinear soft SUSY breaking scalar terms, 
respectively, given by 

(t^A = arg(AM*/2) <Pb = s^rgiBM*/,) , (1.3) 

where M1/2 is the universal gaugino mass. The problem is that from 1-loop diagrams in- 
volving squarks, these phases must be fine-tuned to order 10~^ - 10^'^ to satisfy the limit on 
the neutron electric dipole moment unless all the superpartners are "heavy" , ~ 1 TeV. With 
CP spontaneously broken in a sector independent of SUSY breaking, these phases would be 
naturally zero at first order. 

Attempts to realize the Nelson-Barr mechanism in SUSY models have, however, run 
up against a formidable difficulty: There generically exist potentially large 1-loop contri- 
butions to 6 in these models [|l8l. The dangerous diagrams are shown in Fig.l, where now 
in the supersymmetric case 6 also gets contributions from the argument of the gluino mass 
(Fig.lb): 

9 = e-aTg detMq - 3 argM,, . (1.4) 



As discussed at length in Ref. [|18|, one requires an exceptionally high degree of proportional- 
ity of the soft SUSY breaking trilinear scalar couplings to their associated Yukawa couplings 
as well as degeneracy among the soft squark mass terms for each charge and color sector, 
if these contributions are to be sufficiently suppressed. This is equivalent to the statement 
that when the quark and squark mass matrices are diagonalized by the same set of unitary 
matrices, no phase can appear in the diagrams of Fig.l. The degree of proportionality and 



degeneracy required among the soft SUSY breaking parameters is very difficult to maintain 
due to the effects of renormahzation. 

So, is SCPV doomed to be disfavored as a solution to the strong CP problem in su- 
persymmetric models? We will argue that models with a specifically modified Nelson-Barr 
mechanism together with the recently popular gauge-mediated SUSY breaking (GMSB) sce- 
T9|^3| can overcome the difficulty. The GMSB scenario ensures that the soft masses 
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at the intrinsic SUSY breaking scale Mmess are proportional and degenerate, while renor- 
mahzation effects that violate these conditions are reduced by having to run soft masses and 
couplings from the messenger scale M^ess — 10 — 100 TeV instead of the reduced Planck 
mass ~ 2 X 10^*^ GeV. In the minimal version of such GMSB models [^,^, the A- and 
S-terms are zero at Mmess- This tends to give additional suppression of the dangerous 
contributions to 9. However, large third generation Yukawa couplings can still lead to sig- 
nificant violations of proportionality and degeneracy and a detailed numerical analysis of 
the situation is necessary to determine if the supersymmetric Nelson-Barr type models are 
viable solutions to the strong CP problem. We will answer this in the affirmative with what 
to our knowledge is the only full renormahzation group (RG) analysis of such models in the 
literature. 

The models to be discussed here have CP spontaneously broken at low energies (of order 
a TeV), and a source of weak CP violation distinct from that in the Standard Model, namely 
the exchange of a new U(l) gauge boson, the symmetry of which enforces the Nelson-Barr 
texture. The KM phase is very small. This makes it possible to account for the smallness 
of 9 without making weak CP violation inconsistent with observations. In the SUSY-GUT 
Nelson-Barr models, even with gauge-mediated SUSY breaking, one predicts a too large 9 
if the experimental requirement that 5km ~ C^(l) is imposed. Moreover, one has a richer 
and quite distinct phenomenology. The extra quarks and other fields needed to construct 
the Nelson-Barr texture and break CP will be within the reach of future accelerators. A 
non-supersymmetric version of the type of model is the aspon model [P^-pB| . The situation 



with SUSY incorporated is first discussed in Ref. [^] where its advantage over the generic 



SUSY Nelson-Barr type models is highlighted. Another possible advantage for such low scale 
models would arise if it were somehow possible to imbed the sector responsible for the Nelson- 
Barr texture in the SUSY breaking and messenger sectors. An immediate disadvantage of 
this approach is that by breaking CP at low energies, one introduces a serious domain wall 
problem However, this can be solved via a period of inflation just above the weak scale. 
Indeed, in the SUSY context, some authors have argued that this type of inflation can be 
natural and desirable for other reasons {e.g. as a solution to the cosmological moduli 
problem) . 

This article is organized as follows: In Section 2, the considerations leading to viable 
models for solving the strong and supersymmetric CP problems are discussed. We also note 
some intriguing alternatives worthy of further consideration. We give in Section 3 a detailed 
summary of the dangerous 1-loop contributions to 6 and their dependence on proportionality 
and squark- degeneracy violating mass insertions. The question of the detailed structure of 
the spontaneous CP breaking part of the superpotential is taken up in Section 4. We 
emphasize its importance in determining certain dangerous contributions to 6 and present 
a minimal example we use in further analysis. The renormalization group analysis used to 
estimate the 9 contributions is described in Section 5 and the full numerical results presented. 
Some remarks on related questions of interest are presented in the conclusion. 



II. MODEL-BUILDING CONSIDERATIONS 

We assume full CP symmetry in the visible sector, including the soft SUSY breaking 
terms, down to energies where spontaneous CP breaking occurs. From Eq. (|1.4| ) we see 
that there should be no tree-level phases in the quark mass determinant, nor in the SUSY 
breaking gluino mass. The latter holds by assumption and the former is obtained through 
a Nelson-Barr texture. To obtain the texture, we introduce an extra heavy right-handed 
down quark superfield D together with its mirror D coupling to the ordinary down quarks 
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via the superpotentialQ 

Wd = Y'jQ.m, + ^^,DD + f "Dxa^. , (2.1) 

where the VEVs of the scalar components of Xa contain a relative phase, thus breaking CP. 
The details of the superpotential accomplishing this are postponed to Section 4, and here it 
is sufficient to note that at least two x fields are necessary, since one phase can always be 
absorbed by a field redefinition of the extra quarks. After CP and SU(2)i4/ x U(l)y breaking, 
we have the down sector fermion mass matrix: 



ruq 



(2.2) 



y fio J 

where is the usual 3x3 down sector mass matrix and a is a complex 3-vector with 
components a* = ~^7*"(Xa), and the real parameter x is defined such that a is normalized 
to 1, i.e. a^a= 1. The magnitude of mixing between the ordinary quarks and the extra 
singlet is characterized by x. Clearly the determinant of nig is real and at energies below /i^ 
the low energy effective theory has a KM phase of at most order x. Without some additional 
source for weak CP violation, this must be 0{1) and as we shall see this in turn makes the 
suppression of 1-loop contributions to 6 problematic.^ For this reason, we shall assume that 
CP is broken at relatively low scales with a nonstandard mechanism for weak CP violation. 

The specific form of the mass matrix in Eq. ( ^.21 ) can be enforced by a variety of symme- 
tries, though by the non-renormalization theorems, the Nelson-Barr texture is not upset by 
renormalization of terms in the superpotential. In the aspon scenario, as discussed in Ref. 
||25| , the D, D and Xa can be given charges 1, —1 and 1, respectively, under a new gauged 
U(l) symmetry. The major source of weak CP violation, in the K — K system for instance. 



^ Other type of phenomenological features from an extra vector-like quark singlet have been 
studied by various authors in a different context. See for example Ref. |^0|, and references therein. 

^This is exactly the reason why the small x scenario is discarded in the analysis of Ref. 
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then comes from exchange of the new U(l) gauge boson (aspon) which becomes massive at 
the scale where CP is spontaneously broken. This places an upper bound on the mass scale 
of CP breaking of 0{TeV). More important in our SUSY version, this allows the parameter 
X to be small, e.g. ~ 10~^, which contributes significantly in suppression 9 from loop 
corrections. Note that we will need at least some extra mirror partners for Xa superfields to 
cancel the gauge anomaly introduced by their fermionic components. 

If the Nelson-Barr texture is obtained from a discrete/global symmetry, one must rely 
on superbox diagrams involving gluino and chargino exchange to generate ex'-, a scenario 
recently re-analyzed in the context of the minimal supersymmetric Standard Model (MSSM) 
||3T| . However, in present setting, this proves very difficult to do. This is unfortunate since 
it is easier to construct unifiable models in this latter case.0 

When SUSY is broken we generate nonzero gaugino masses and soft scalar bilinear and 
trilinear couplings]^ including the following terms relevant for the down-sector analysis [a 
complete description is given in Eqs.(^,^)]: 

Koft = Q^^lQ + dmli + b^mlb + bmlD" + Xa^lahXl 

+ dhdQHd + Ihi^bxa + Bu^luDb + h.c. (2.3) 



The general form of the down squark mass matrix can be written as. 



(2.4) 



To construct a unifiable model in the discrete symmetry case, we need to add a pair of heavy 
lepton doublets coupling in a superpotential analogous to Eq.(^j|). The same thing can be done 
to make the aspon model GUT-compatible, but unifying the extra U(l) with the other gauge 
interactions is not possible. 

^Note that the trilinear couplings, and /i!^", are also commonly written as products of the A- 
parameters and the corresponding Yukawa couplings, e.g. elements of the hd matrix correspond to 
■^d'^d sum). 
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where 



Ml, 



/ ~2 t 



x^la^ 



hdVd + mdidH tan(3 Mjb ^ 



In the expression for A4%l '^^ have used 



(2.5) 



(2.6) 



(2.7) 



Mlh'^hl%Xa)-r{F^:) , 



(2.8) 



where b is normahzed to 1 and F^^ is the F-term for the Xa field, which depends on the 
specific form of the soft SUSY breaking mass terms related to the spontaneous CP breaking 
part of the superpotential. The form of these squark mass matrices will be critical in 
the calculation of 9 in the next section. In particular, the {Fy,J^s bear complex phases 
independent of those in the (Xa)'s in a generic setting, and hence constitute a major source 
of trouble. 

We have implicitly assumed in the above that the sectors responsible for the Nelson- 
Barr texture and CP breaking are disjoint from those involved in the intrinsic breaking of 
supersymmetry. Since the successful example model we focus on has gauge-mediated SUSY 
breaking at a relatively low scale Mmessi it is a priori possible that all or part of the extra 
field content and symmetries required for the Nelson-Barr mechanism is contained in the 
SUSY breaking hidden and messenger sectors. Although this is an intriguing possibility, we 
have not been able to construct viable models of this type thus far. Actually, there appears 
to be an intrinsic incompatibility between the role a field takes in CP violation and the one 
it takes in SUSY breaking, as far as constraining 9 is concerned. Recall that in GMSB, the 
scalar particles get their soft SUSY breaking masses from the gauge interactions they share 
with the messenger sector particles which see SUSY breaking directly. The squark masses in 
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each sector, for instance, would then be degenerate, as the process is flavor bhnd. The extra 
singlet D introduced here may easily upset the situation. Naively, the best strategy is to 
make the GMSB also blind to the aspon U(l). We will see below that this happens to have 
a even more important merit — it guarantees the suppression of the very dangerous {F^^)^s. 
In other words, hiding the CP-breaking sector from SUSY breaking helps to suppress the 
SUSY loop contributions to 6. This is the less ambitious strategy we have taken in the 
model analyzed in detail below. 



III. CALCULATION OF 1-LOOP 6 CONSTRAINTS 



Here we review and extend the work in Refs. |18,27| to compute the 1-loop contributions 



to 9 of Figs. la and lb, using the mass-insertion approximation These results are 

generally valid for any type of low energy SCPV model, with or without the U(l)^. The 
scale, fijj, and hence the characteristic scale of the SCPV, here is chosen near or below 
•fhsq, the average squark mass. This is more or less dictated by the aspon scenario of weak 
CP [P5[. The D and D superflelds are handled on the same footing as the other quark 



superflelds. The analysis is basically the same as that given in Ref. ||2^ except here we pay 
full attention to the explicit phase factors and family indices, and also treat the M| term, as 
given by Eg. (|2.8|) , in full detail. These turn out to be very important in understanding how 
the scenario can provide a feasible solution. In taking the large fir> limit, which corresponds 



to situation discussed in Ref. [|T^, one has to be careful in handling the loop momentum 



integrals properly. The latter are however not explicitly given in this paper, though they 
are included in our numerical computations. 
The 1-loop contribution is given by 



56 = Im Tim-^Smp + SM'^SM, 
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where m^^ runs over the four eigenvalues of the quark mass matrix [cf. Eq.(p^.2|)] and Ml 
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over the eight eigenvalues for the squarks [cf. Eq.( |2.4|) ], all in the down-sector; Z^'^ is the 
unitary rotation that diagonalizes the squark mass matrix in the quark mass eigenstate basis. 
This full formula, while it can be used in the numerical calculations once all the quantities 
involved are known, hides its physics content behind the Z-matrix elements. In the limit 
of exact degeneracy and proportionality, the latter is just the identity matrix and Q is zero. 
Otherwise, the mass-insertion approximation, as discussed below, is more illustrative. 

We first assume an approximate degeneracy and that the diagonal blocks in A^^ dominate 
over the off-diagonal block A^^^ and write 

rrig = frig x 1 + 5m j , ^n^j = frig + om^, , 

= ""^d X 1 + om^ , = ffi^ + om^ . (3.2) 

The squarks are then treated as scalars of masses m| and with the 5fn^ and 
treated as admissible mass-insertions is the loop-diagrams Figs, la and lb. Explicit forms of 
the matrices needed to diagonalize mp are useful. Expressions up to order are available 
in the literature . To parametrize the effect of proportionality violation among the three 
families, we write 

hd = AaYa + 5Ad . (3.3) 

The situation for the related parameter in the d-D mixings is more complicated. Recall that 
MlU = h}^{Xa) - Y'^iFxa) [Eq-dl]!)]. It has been emphasized in Ref. p3 that the F-terms 
being small is paramount to the success of any model of the Nelson-Barr type. These terms 
are dangerous because in general one has no reason to expect these F-terms to obey even 
an approximate proportionality (to the x/ioO* terms). On the contrary, contributions of the 
other part to 6 can be interpreted as a proportionality violation among the 7*'^'s by writing 

= A^Y'" + ; (3.4) 

the term proportional to A^ does not contribute. We further introduce the simplified nota- 
tion: 
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5A^d = /i!; (xa) - A^a' , (3.5) 

where complex vector c is normalized to 1. Hence, we have 

= A^{xfina') + 5A^{xfXnc') - (F^J . (3.6) 

In terms of the above notation, the list of major contributions to 6 is given in Tables la 
and lb. The 6 contributions involving M| are complicated. To make it easier to see the 
effects of the different parts, we list some of those terms in tables before and after the above 
mentioned splitting. For example, entry 1 in the Table la is split into two parts: the first 
part is a proportionality violation effect involving 6A^ and Im(a*c*) (both are suppressed in 
our model), the second is the F-term contribution (7*" (F^^)), where the relevant complex 
phase is taken to be 0{1). One other notable feature among the 6 contributions is the 
combination M|6* — xfioBod^ as shown in entry 9 of Table lb. When the M|6* term is 
split as above, the second term actually can be combined with the first term in Eg. ( p. 61) 
to give xfio{A^ — Bo)a} 1 which can be interpreted as a proportionality violation among the 
corresponding trilinear and bilinear terms. The other parts involve 5A^ and 7*" {F^a)i ^ 
explicitly shown in the table. All other entries with a M| can be split and interpreted in the 
same way. We will see in the final result that the F-term contribution is the most dangerous. 



IV. THE SPONTANEOUS CP VIOLATION SECTOR 

Spontaneous breaking of the U(l)^ symmetry is the only source of CP violation in our 
model. This CP violation effect feeds directly into the x/i^ and M| terms in the quark and 
squark mass matrices, with complex phase vectors a and b, respectively. To implement the 
mechanism, we need a sector of U(l)^-charged SM singlet superfield with a superpotential 
that not only gives rise to the complex (Xa)'s, but also gives us a good control on the 
dangerous (F^^)'s. Soft SUSY breaking terms should also be taken into consideration, when 
determining the true scalar potential. The F-terms, of course, characterize SUSY breaking. 
We consider the scenario in which the messengers communicating SUSY breaking to the 
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visible sector are U(l)yi-blind, i.e. they do not carry any U(1)a charges; furthermore, they 
are not directly coupled to the CP-breaking sector. The superfields of the latter are then 
hidden from SUSY breaking. 

We have to consider at least five superfields, two x's of conjugate U(1)a charges to the 
Xa's and a singlet K, in order to have both gauge anomaly cancellation and a possible CP 
violating vacuum solution |l27| , |3^ . We consider the superpotentialQ 



= Xal^l'Xb + ^XaX'^'Xb + K^' + f^.^' ■ (4.1) 

The five F-flat conditions yield four independent equations, which, together with the D- 
fiat condition, give a unique vacuum solution. The solution is CP violating for most of 
the parameter space. Hence, neglecting the soft SUSY breaking terms, we have a SUSY 
preserving vacuum that breaks CP. 

The GMSB scenario we considered allows the unwanted soft SUSY breaking terms of 
the sector to be zero at M^ess- They are, however, generated through RG evolution, as 
discussed in the next section. With their nonvanishing values taken into consideration, the 
scalar potential is then given by 

Vx = Dl + F^Fl + F^Fl + F^F: + (4.2) 

where 

+ xiKabXb + XaKabXl + ^^K^ ■ (4-3) 

Solving for the potential minimum to determine the (F^^) values is not tractable, as 
and Vsy involves a large number of parameters which are not otherwise constrained, apart 



^In Ref. 1 27], a without the /i^-terms is suggested. While that could have a CP violating 
vacuum with SUSY preserved, the situation is not as general and natural as the one considered 
here, and would have to rely on a linear ^!; term to fix the symmetry breaking scale. 
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from yielding a CP violating solution. However, one can easily obtain a reasonable order of 
magnitude estimate of the shifts in the (F^^)'s as a result of including the small Vs^ terms. 
For example, the equation 

— ^ = 2D^-^ - (/if + ^\^')F^, + F,-^ + B^'x, + hfxb^ + xlrhl,^ = (4.4) 

CXci ^Xa OXa 

suggests that {F^ (here we drop all indices and phases) is given by the magnitude of 

or ?fi| or h\ {x) /A . (4.5) 

An alternative way to estimate the {F^^)^s is given by the SUSY breaking diagrams shown in 
Fig. 2. Here Figs. 2a and 2b give exactly the same results as the first two terms listed above. 
Figure 2c, however, gives the estimate (F^) ~ hx {x) A/lGvr^. For perturbative values of the 
A coupling, this is of course smaller than the third estimate in Eq.( [4.5|) , hence we neglect 
it. A similar diagram, Fig.2d, also suggests a contribution ~ /i^ (%) 7/167r^, though the 
7 dependence of the {F^)^s is implicitly incorporated into the generation of the Vs^ terms 
through RG running. We will use all these in our numerical estimates to determine whether 
the F-term is sufficiently small that its contributions to 6, listed in Tables 1 and 2, are under 
control. Finally, we emphasize again that the complex phases in the (F-^^)'s are not related 
to those of the (Xa)'s directly. 

V. RENORMALIZATION GROUP ANALYSIS 

As pointed out in the introduction, we need a full theory for the soft SUSY breaking 
parameters to see if the 6 constraints can be satisfied and the GMSB scenario may provide the 
only viable possibility. In particular, we use here only the minimal version of such a theory 
pO| . This version has a few special merits: it provides practically a one-parameter model of 
soft SUSY breaking, radiative breaking of electroweak symmetry is naturally implemented 
and, within the MSSM framework, it has been studied with extensive renormalization group 
analysis and shown to be compatible with all known experimental constraints p2|. From our 
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perspective of solving the strong CP problem by augmenting the U(l)^ sector, it actually 
represents a relatively demanding setting among GMSB models, where a large tan/? allows 
all the Yukawa couplings of the third family to have substantial effects on the RG-runnings. 
A smaller tan/5 in general would only make it easier to satisfy the 9 constraints. 

We will refrain from elaborating extensively on the details of the GMSB model or the 
RG-analysis itself. For more specific details on finding the correct electroweak symmetry 
breaking vacuum and meeting other experimental constraints in the minimal GMSB model, 
readers are referred to Ref. [^. Our interest here is in adapting the machinery to our 
extended model at a level of sufficient sophistication to calculate 9 to 1-loop and establish 
our solution to the strong CP problem. 

We use 1-loop renormalization group equations (RGB's) with naive step thresholds be- 
tween Mz = MsusY and Mmess- The RG-improved tree level Higgs potential is considered 
in finding the electroweak-symmetry breaking solution. The RGB's for the extra content of 
the model are also implemented at the 1-loop level; relevant formulae are in Appendix A. 
The U(l)^, and hence CP symmetry, breaking is imposed by hand. The idea is to study 
the general situation independent of the details of the SCPV sector, as the latter is to a 
certain extent more flexible and less constrained. It is important to note that the extra 
superfield content in the model is partially decoupled from the MSSM part, with the only 
direct coupling being gauge couplings of D and Z), and the small Yukawa couplings 7*". 
The computation concerning the SCPV sector, as well as the RG analysis can certainly be 
made more sophisticated, however, we consider our treatment sufficient for our purpose. In 
the sample analysis for which numerical results are presented in this paper in detail (Ap- 
pendix B and the last column of Table 1), the values of the various 7*" Yukawa couplings are 
generated randomly in the range 0.005 — 0.01. The latter is chosen to target an x-value of 
around 0.01. The value of yU^, is fixed at 500 GeV; Mmess ~ A at 50 TeV. For the soft SUSY 
breaking parameters from GMSB, all A- and B- terms are taken to be zero at M^ness = 
The scalar soft masses from GMSB are given by 
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m\X) = - {c,al{X) + C,al{X) + f{y) , (5.1) 

where C3 = 4/3,0 for triplets and singlets of SU{3)c, C2 = 3/4,0 for doublets and singlets 
of SU{2)l] F = Q — T3 is the hypercharge. The function /(y), derived in Ref. [^, is 
simply set to 1. Note that the above formula is independent of the U(l)^ charge; a SM 
singlet with or without U(l)^ charge, such as the Xa and K scalars, has no initial soft mass. 
Gaugino masses are likewise given by the MSSM formula, omitted here. The new U(l)^ 
gaugino (aspino) has no tree-level SUSY breaking mass. The gauge coupling qa is taken to 
be around gem- The SUSY-breaking aspino mass Ma then remains vanishingly small even 
after finite loop effects and RG evolution are taken into account. 

After the electroweak symmetry breaking solution is obtained, various 9 contributions 
are calculated through the mass-insertion approximation to order x^. To impose the U(l)^ 
symmetry breaking, we set, in the sample analysis, | (xi) P + | (X2) P — and choose 
random values for the VEVs and their complex phases within the constraint. Effects of 
higher order in x are checked to be insignificant. 

Values of parameters in Wy^ are needed for the RG-runnings of particularly the SCPV 
sector soft SUSY breaking parameters, discussed in the previous section. To simphfy the 
situation, we input all these mass parameters as /i^ and all dimensionless couplings as 
random numbers in the range 0.1 — 0.8, for the sample calculation. This oversimplification 
certainly begs the question of consistency of the vacuum solution for this sector, or the whole 
model. However, in the small x domain of interest, the influence of the extra ingredients 
on the values of the other MSSM parameters is insignificant, as to be expected. The only 
practical effect of those parameters is in the RG evolution of the related soft terms which 
we needed to estimate the {F^)^s. We have checked, for instance, that the particular input 
values used in the sample run reported here does lead to generic magnitudes of the latter. 

Appendix B contains a collection of some of the numerical results, while those for the 
6 contributions, without the {F^)^s are listed in the last column of Table 1. Estimates of 
the (F^), following the discussion in the previous section, and their contribution to 6 are 
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given in Table 2 (second column). The latter can be easily checked using the (F-^^) value 
and the listing of M|-terms in Table 1. We also Hst in Table 2 results from a number of 
different runs with different values of the 7's (reflected by the x- value obtained) and A's. The 
former, which can have a significant effect on the various MSSM parameters, are restricted 
by ~ 10~^ - 10~^ from the weak CP considerations. Our results indicate that only a 
relatively large value of x can upset the strong CP solution, by first driving (F^) too large 
(see column 4 of Table 2). One should be cautious in using this result quantitatively, as our 
(F^) estimates are meant to be conservative upper bounds. However, the result is certainly 
illustrative of the importance of the (F^) in estimating 9. With x restricted to the workable 
range, the basic features of the RGE solutions are quite stable. This is true even with a 
relatively large variation of the A's, as illustrated by column 5 and 6 of Table 2. Note that 
though the hx (%) /A term may have an explicit dependence on A, its numerical value does 
not have a large variation with A as one might naively expect. This is, like the approximate 
proportionality of a general A-term, a natural result of the RG equations. 

All in all, the F-term contributions to 6 dominate, and the overall 6 value is comfortably 
within the required bound for the major region of the parameter space of our model under 
consideration, hence solving the strong CP problem. 

VI. CONCLUSIONS 

To recapitulate, we discussed a complete spontaneous CP violation model with gauge- 
mediated supersymmetry breaking and why this type of model is particularly favored over 
a generic supersymmetric SCPV model in solving the strong CP problem. Results from 
numerical RGE studies are used to exphcitly estabhsh the feasibility of the approach. The 
treatment of parameters in the SCPV sector is admittedly oversimplified. The correlations 
between x and Ho, and between the various mass parameters at the /lo scale and the values 
of the various A coupling, for instance, are neglected. However, it is easy to see from our 
discussion that such details are not going to change the essential features of our results. 



18 



though they would determine exphcitly the specific "large-x" region of the parameter space 
that could be ruled out. The model has a rich spectrum of new particles at the fij^ or SCPV 
scale. Until such experimental data become available, a detailed study of the parameter 
space may not be feasible. 

While the weak CP aspects of this model have been analyzed in the non-supersymmetric 



setting as in Ref. pSf, SUSY particles could lead to new contributions through super-box 
and penguin diagrams. These contributions are in general subdominant, as are the Standard 
Model box diagrams. 

Our model predicts a measureable neutron EDM, which could be close to the present 
experimental bound for x > .01. The model also has an extra pair of vector-like quark 
superfields, a new gauge boson, and a number of neutral fermions and scalars with no direct 
couplings to the Standard Model gauge bosons, all with masses around the TeV scale. This 
scale is dictated by the weak CP phenomenology. Hence it offers a rich spectrum of new 
particles to be discovered at the CERN Large Hadron Collider. Because the 3x3 Cabbibo- 
Kobayashi-Maskawa matrix is approximately real, the model also predicts a flat unitarity 
triangle and the absence of substantial CP violation in the B system at future B factories 
2^. Moreover, there will be a lack of any substantial CP violating effects in the up-quark 



sector. 



ACKNOWLEDGMENTS 

The authors want to thank P.H. Frampton for being a constant source of support and 
encouragement. O.K. is in debt to colleagues in Rochester, where the major part of this 
manuscript was finished. Special thanks go to M. Bisset for reading the manuscript. B.W. 
thanks P. Sikivie for discussions. The authors were supported in part by the U.S. Depart- 
ment of Energy under Grant DE-FG05-85ER-40219, Task B; O.K. was also supported in 
part by the U.S. Department of Energy under Grant DE-FG02-91ER-40685. 



19 



APPENDIX A: RENORMALIZATION GROUP EQUATIONS 



Here we collect the modified one loop MSSM renormalization group equations (RGE's) 
to account for the extra vector-like chiral superfields and the new Yukawa couplings in our 
model. In many cases we give only the extra contributions and refer the interested reader 



to Ref. I^Sf with whom we share conventions. The complete two-loop renormalization group 



equations for a softly broken supersymmetric theory can be found in Refs. p5| , |36 
The complete superpotential can be written as 

+ Xa/i?X6 + ^XaX'^'Xb + AkK3 + (Al) 

where family indices are implicit except in the new Yukawa coupling and a, 6 = 1,2. The 
soft super symmetry breaking Lagrangian can be written as £soft = -^^ft^ + where 



= uKQHu + dhdQHd + eKLHd + B^fi^H^Hd + h.c. 

+ Q^rhlQ + Dfh\L + uffilv} + dffili + erhle^ (A2) 
+ ml^HlH^ + mlHlH, 



and 

- Ct^l^ = d,h\^DXa + B.fi.DD + X + Xahfxb^ + + B^^^ + h.c. 

+ bmlb^ + b^mlb + xlrhlabXb + Xa^^labxl + ^^K^ ■ (A3) 

Note that B^ and B^ are defined in a different way from i?^ and Bu] the former have 
dimension (mass)^ and are analogs of Bj^fi^ and Bu^g. Also, 5^ is a 2 x 2 matrix. Finally 
we have supersymmetry breaking gaugino masses Ma{a = 1,2,3) for the MSSM and a 
possible gaugino mass M4 under U(l)^. 

We give the MSSM one loop RGE's for these interactions below. The gauge couplings 
are computed to two loops and are given by 

V / \b=l x=u,d,e,'y J 
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where = (f + f + |7V^, 1 + N^, -3 + No), 



( 



5(2) 



199 I 8 AT- I 9 AT- 27 I 9 a/- 88 , 32 at- 



5 ' 5 



and 



C 



25 + 7A^L 
9 



26 14 18 4 

5 5 5 5 

6 6 2 
4 4 2 
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14 + ^Nd 



(A5) 



(A6) 



In the above we have allowed for the possibihty of heavy vector-hke pairs of charge | 
color triplets [SU(2) singlets] and such pairs of hypercharge —1 SU(2) doublets (color 
singlets). These include both the extra mirror pairs D + D and L + L which can interact 
directly with MSSM fields, but also extra mirror pairs originating in the messenger sector at 
higher scales. We have for simphcity omitted the possible Yukawa interactions of the extra 
mirror lepton doublets and the two loop U(l)^ contributions. In the actual computations, 
we use No — 1, Nl = 0. At one loop we also have 

^dgA 



167r" 



dt 



(A7) 



when the U(l)^ is present. 

Using the above definitions, the two loop gaugino mass equations are 



dt 167r2 " (167r2)2 



(A8) 



At one loop we have 



167r^ 



dt 



20glMA . 



(A9) 



The running down quark Yukawa matrix is modified by the presence of the D-d, couplings 
7*" which are treated as 3 x 2 matrices below. We have for the up and down one loop Yukawas: 
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igtt^^ = Y, (Tr (srir, + ntyj + srir, + Kjy; - ^^3' - s^?^' - ^gl) + 77^>^. , (aio) 



dt 
,dY,, 



levr^^ = Y. (31V (yjy„) + sy^y^ + yjy, - - ^gl - ^^gi) , 

and for 7: 

IGtt^^ = 7 (tt (7^7) + 27l"7 + AU 



(All) 



(A12) 



where for 7 we have included the effect of a possible extra U(l)^ as described in the text. 
The running supersymmetric /x^, parameter is given by 
id^D / + \ 16 r, 4 



The equation for the corresponding soft mass parameter is 



,dB 



IGtt^^ = 2Tr (7^/.,) + ^-^glMs + :^^,^M, + S^^M^ 



dt " ^' "'^ ■ 3 

Similar parameters for the SCPV sector have 
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IGtt'^ = (aU + 37^7) + AA^S^ + A [Tr (a^sJ + X^B^ 



dt 



and 



IGtt 



167r^^ = /.. [A^, + 2IV (At a) 



= [A^ + 2Tr (aU)] + A, [A,S, + 2Tr (A^eJ] 
+2/X, [Ah/i. + 2Tr (At/z;,' 



There are also RGE's for the extra Yukawa couphngs: 
ydX 



16n^^ = A (iV (AtA) + ^A^ + 2AtA + 3^7 - ^gl) , 



levr^^ = A. [A^, + 2Tr (AtA) 



dt 



(A13) 



(A14) 



(A15) 



(A16) 



(A17) 



(A18) 



(A19) 
(A20) 



22 



Note that fi^ ,X and hx are all 2 x 2 matrices. 

The relevant soft supersymmetry breaking trilinear coupling RGE's are given by 

levr^^ = h, (Xr {3YjY, + Y^) + 5YJy, + Y^Y^ - fgj - Sgl - ^,?) + 

+ yJtt (QYjh, + 2y^k) + ^ylK + '^ylK 



+ y (?32M3 + QglM2 + ^glMi) + 2h,^^Yd , (A21) 
16vr^^ = K (sTr (fJfJ + 5^^^^ + Y),Y, - ^-^gl - 3^,^ - ^^gl 
+ i;(6Tr (y>J + 4F>„ + 2Fi/i, 

+ y + QglM^ + , (A22) 

levr^^ = ^ (Tr (7^7) + 37^7 + AtA - ^^3^ - - 4^i) + Y.Y^h, 
+ 7 (^2Tr (7+/^^) + 37^/i-, + At/iA + ^glM^ + + 
+ 2/i,yi7, (A23) 



/lA f Tr (AtA) + 37^7 + SA+A + ^A^, - Ag\ 



+ A (2Tr (A^/ia) + 67^/1^ + 3At/iA + A,/i, + SglMA) , (A24) 
IGtt^^ = (^^Xl + 3Tr (A+a)) + GA^Tr (X^h^) . (A25) 

Finally we give the modifications of the MSSM soft hermitian quadratic mass parameters 
(see Ref. for a complete description). All of the MSSM RGE's have the following change 
in a D-term contribution: the factor S defined in Eq.(4.27) of |^ is now given by 

^ = - -mjj + m% + TT (ml -ml- 2m? + m| + mfj . (A26) 

Besides this modification, the equation for m| has the only nontrivial change due to the 
coupling 7: 

16^2^- ^ ^2m| + 4m^J Y,Y} + AY.mlYj + 2Y,Yjm} + 77^m| + rh^fl^ 
+ 2m^77^ + 27mj7^ + Ah^hl^ + 2h^h\ 

-fglmf-^^glm' + lglS, (A27) 
where is a 2 x 2 matrix. The equations for the new soft squark masses are 
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tdrh 



dt 



: 2Tr (^m^7^7 + 77/1^7''" + 7^m|7 
s 



167r^ 



^gl \Msf - ^gf \M,f - %g\ \Ma\^ - '^-glS + 2g\Sj, 



D 



dt 

dm? 

A_ 

dt 



32 2,,. 



10 ; 

2 



^dml 

dt 
,dml 



: 37 '7m^ + 3m^7'7 + Gj^^rrij-, + 67 'mj 

-2 



m|7 + 6/i;^/i^ 



(A28) 
(A29) 

(A30) 



■ '2g\SA , 



dt 



i\ + mjAU + 2AUm^ + 2Atm|A + 2^/1;^ -%g\\MAf - 
m\\\^ + AA^m| + 2AA^m^ + 2AmjA^ + 2h^h\ - Sg\ \Ma\^ + 2g\SA , (A31) 
?^ + 2TrAU) ml + 2Tr (AmjA^) + 2Tr (A%|a) 



= (3A^ + 

+ li. 



(A32) 



where 5yi = 3m|, — 3m|| + Tr(m| — m^). 



APPENDIX B: SOME NUMERICAL RESULTS 

We collect here some of the numerical results in our sample calculation. Recall that we 
use Mmess = 50 TeV, fx^ = fi^ = 500 GeV, and random values of 7's and A's in the ranges 
.005 — .01 and .1 — .8 respectively. As inputs we also used the values Mt — 175 GeV and 
as = 0.12. 

The electroweak-symmetry breaking solution is obtained with 



tan/3 = 43.18 , i^h = -370.5 GeV , Bh = 3.938 GeV . 



(Bl) 



In the following, we concentrate on the down-sector as it is the only one of relevance to 
the understanding the 9 value. Average values of the squared left- and right-handed squark 
masses are 5.337 x 10^ GeV^ and 5.031 x 10^ GeV^, respectively, and is -243.5 GeV. 
The lack of proportionahty of the Ad-terms is given by 

^ 0. 2.12 X 10-^ 8.45 X 10"^ ^ 

1.18 X lO-'^ 6.39 X 10-^ 4.02 x 10"^ , (B2) 
-2.37 x 10"^ 1.73 x 10-6 1.63 x lO'^ 



5A/m 



sq 



\ 



I 
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while degeneracy violations are given by 



and 



\ 



0. 2.54 X 10^5 -6.13 X 10"^ 

2.54 X 10-5 -2.53 x 10"^ 4.42 x 10"^ 

-6.13 X 10^^ 4.42 X 10^3 -0.172 ^ 

0. -8.22 X 10^^ -1.12 X 10^5 ^ 

-8.22 X 10-6 -1.59 X 10-^ -1.17 x 10"^ 

-1.12 X 10-5 -1.17 X 10-5 -0.135 



where m^g is the average squark mass. In addition, we have 



Srhl/rhl = 1.65 x lO'^ , 



6mj)/mig = -5.64 x 10" 



(B 



(B 



(B 
(B 



The proportional part of the A-terms for the 7 Yukawas is very close to the fi/j value, giv 
by 



A^r^ Bd = -248.1 GeV 



(B 



with a difference of only 4.27 x 10 GeV. Some other quantities of interest are 



MJ^ M,) = 478.1 GeV , x = 1.213 x 10"^ 



and, as defined by Eg . (p75| 



SA^ = 4.782 X 10-^GeV , Im(atc) = 5.531 x 10' 



fB 
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Table Captions. 

Table 1: Analysis of the major 1-loop 9 contributions and numerical results from the 
sample run. Table la contains contributions from gluino mass corrections; Table lb contains 
those from quark mass corrections. Entry 1 of Table la and entries 9, and 10 of Table lb 
are shown together with explicit splittings of M| according to Eq.(3.6) below the first lines. 
Numerical results given in the last column do not include the {F-^ term contributions, 
but otherwise are complete, i.e. they include all other numerical factors from color indices 
summation, momentum loop integrals, and full summation over family indices (i, j, k) so that 
the full 9 value without the F-term contributions, apart from some unlisted subdominating 
terms, is given by the sum of all the entries. 

Table 2: Estimates of the (F^) term and its contribution to ^, for our sample run and 
a few runs with different 7 and A inputs {^u and //^'s are all set at 500 GeV, M^ess at 50 
TeV). Note that the entries B^, m|, hx{x) and /i-y (x) 7/1671^ are our (F^) estimates, 
as discussed; all these are quantities of dimension (mass)^ in units of GeV^ (not shown 
explicitly). The (F^) estimates and its contributions to 9 are meant to be upper bounds. 
Overall 9 contributions from gluino and quark mass corrections without the F-term are also 
listed. 

Figure Captions. 

Fig.l : 1-loop mass-correction diagrams leading to 9 contributions, (a) 1-loop quark 
mass; (b) 1-loop gluino mass. 

Fig. 2 : Diagrams giving estimates of (F^) magnitudes. Note that (%), (x), /x^ and all 
propagator masses in the diagrams can be taken as around the same scale, namely /lo] a 
SUSY breaking vertex or mass insertion is required in each shown. 
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FIG. 1. 1-loop mass-correction diagrams leading to 6 contributions, (a) 1-loop quark mass; (b) 
1-loop gluino mass. 
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FIG. 2. Diagrammatic estimates of (F-^) magnitudes. Note that (xa) and all the mass insertions 
in the diagrams are around the same scale, namely /x^,; a SUSY breaking vertex or mass insertion 
is required in each shown. 
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Table 1: Analysis of the major 1-loop 9 contributions and numerical results from the sample run. 
Table la contains contributions from gluino mass corrections; Table lb contains those from quark 
mass corrections. Entry 1 of Table la and entries 9, and 10 of Table lb are shown together with 
explicit splittings of M| according to Eq.(3.6) below the first lines. Numerical results given in 
the last column do not include the (F^) term contributions, but otherwise are complete, i.e. they 
include all other numerical factors from color indices summation, momentum loop integrals, and 
full summation over family indices {i, j, k) so that the full 9 value without the F-term contributions, 
apart from some unlisted subdominating terms, is given by the sum of all the entries. 



Table la: Gluino mass correction contributions 



No. 








factors of 9 


contribution 








(1) 




IJ-D 

rhsq 




rhsq 

Mq 

rhsq 

Mq 






M| 

M-y 

rhsq 


Im(a*6') 
Im(a*c') 


-5.17 X 10-15 






rhsq 




rhsq 






"^iq 


(cf. table 2) 


(2) 




Mo 


rh\m^ 


rhsq 






Ml 


Im(a*60 


-1.53 X 10-1^ 


rhsq 




Mq 








(3) 


as J2 




"^'iq 


r'risq 

Mq 




rhsq 




Im(a*a*) 


1.53 X 10-1^ 


(4) 


as ™2 


■^q 


fhAm] 


rhsq 


a 

■fn'iq 






Im(a|a') 


1.09 X 10"^^ 


(5) 




rhsq 




rhsq 

Mg 


miq 




miq 


Im(a^6') 


1.27 X 10-1^ 


(6) 


as J2 
4tt-^ 






rhsq 


a 






Im(a*a') 


2.39 X lO-i'^ 




rniq 


Mq 


rniq 


rhsq 




(7) 


as ™2 
47r'^ 


■^q 


vami 
■f^iq 


rhsq 

Mg 




rhlq rhsq 




Im(a^a*) 


-5.99 X 10"^^ 


(8) 


as 
4n 




Vdrrii 
rniq 


rhsq 

Mg 




rhsq 


{Ml? 

"^tq 


Im(6*60 


8.43 X 10-2° 


(9) 


4tt 


rrisq 


VdrnifriA 


rhsq 

Mq 




rhsq 


Mi 


Im(a*a*) 


1.14 X 10-1^ 


(10) 


as rp 

4tt-^ 


rhsq 




rhsq 

Mg 






Mi 

■^q 


Im(a|6*) 


5.26 X 10-22 
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Table lb: Quark mass correction contributions 



No. 








[■< 




( )f 6 c< )nl I'i 


")uli()ii 






(1) 




rhsq 




M„ 

9 

•rhsq 








5 

rh'iq 


Im(a*6*) 


-1.52 X 10-15 


(2) 


Ois „ 


fJ'D 
fhsq 


iiiaBd 
rhsq 


Mg 
Tflsq 








Mi 

sq 


Im(a|6*) 


3.45 X 10-15 


(3) 


Ole Jl 




Vdtmi 
r^sq 


rhsq 






fhsq 




lm{a*a') 


3.10 X lO-i'^ 


(4) 






rrij 


Mg 


Sfii^^^ 

d 






Ml 


lm{a*bj) 


2.02 X 10-1^ 




rrij 


rhsq 


fh'iq 






fhiq 


(5) 


47r'^ 




'"d 


Mg 


5fh?^^ 

d 








Im(a|a-^) 


2.00 X 10-15 


sq 


rrik 


rhsq 


''~"'in 




fhsq 






as T.2 


sq 


Vd 

m-i 
J 


Mg 

rhsq 


ml 

rh'i„ 

sq 




SAV_ 

fhsq 




un.[aj^ a-' ) 


Q fin \/ 1 n— 18 
o.UU X lU 


(7) 


as J.2 
An 


sq 


rriifhA 
nif^rhsq 


Mg 

rhsq 


a 

rni^ 

sq 


fh}^ 

'""sq 






Im{a*a^) 


-7.49 X 10-13 


(8) 


as ^2 


m'iq 


rrij niA 
rrnmsq 


Mg 
msq 


m'iq 


rhiq 








2.26 X 10-13 


(9) 




Mo 
rhsq 


mi 
mi 


Mg 

rhsq 




'"'sq 




Im(a*( 


^6^-XMo|gaO) 


-5.74 X 10-15 




as rf.2 


^ u 


m-i 

3 


Mo 




s^T 




A -Bd 


Im(a*a*) 


— 




rh'iq 


mi 


rhsq 




rn^ 

sq 




msq 








as J2 




rjh. 
mi 


Mg 

rhsq 




rri^ 

""sq 




&Aj_ 

fhsq 


Im(a*c*) 




— 


^X 


rhsq 


mi 


M„ 
rhsq 










rh'iq 


(cf. table 2) 


(10) 


as 
4n 




'^d 

mi 


M„ 
g 

rhsq 






rrtsq 


rhiq 


Im(6*6*) 


1.78 X 10-16 


- 


as '2 
47r'^ 


Ml, 




Mg 










Im(a^a*) 






rh'iq 


rni-miq 


fhsq 






msq 






- 


as J2 
47r-^ 


^<■D 

■^q 


M. 

mi 


Mg 
fhsq 






fhsq 








- 


as J2 






Mg 








5Aj 


Im{c*a' + a*/) 








rriinisq 


rhsq 






msq 


fhsq 






as rp 






Mg 






SA^ 






(cf. table 2) 




rhsq 


miihsq 


msq 






msq 




rhiq 




4tt-^ 


MO 
rhsq 


Vd. 

mi 


Mg 
fhsq 






5A^ 

fhsq 


fhsq 


rn'iq 


(cf. table 2) 




as 
4tt 




M. 

mi 


Mg 
fhsq 






msq 


7 


miq miq 


(cf. table 2) 


(11) 


as 

4-K-^ 


rhsq 


VdBo 
mirhsq 


Mg 
fhsq 






fhsq 


miq 


Im(6*a*) 


3.08 X 10-21 
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Table 2: Estimates of the (Fy^) term and its contribution to 9, for our sample run and a 
few runs with different 7 and A inputs {fio and /i^'s are all set at 500 GeV, M^ess at 50 
TeV). Note that the entries B^, m|, hx{x) /A, and /i-y (x) 7/167r^ are our (F^) estimates, 
as discussed; all these are quantities of dimension (mass)^ in units of GeV^ (not shown 
explicitly). The (F^) estimates and its contributions to 9 are meant to be upper bounds. 
Overall 9 contributions from gluino and quark mass corrections without the F-term are also 
hsted. 



No. 


1 


2 (sample) 


3 


4 


5 


6 


X 


.0081 


.012 


.02 


.068 


.0086 


.0077 


A 


.23 - .44 


.18- .68 


.53 - .77 


.27- .61 


.0038 - .0094 


.75 - 1.3 




1.5 


3.4 


9.5 


101 


1.9 


1.5 


ml 


3.9 


3.5 


.24 


18 


4.7 


2.5 


hx (x) /A 


2.3 


1.0 


11 


7.9 


4.6 


1.8 


167r2 


.033 


.054 


.16 


1.9 


.036 


.032 


(F^) estimate 


10 


10 


25 


130 


10 


10 


' J 2^" estimate 


10-^ 


10-^ 


10-^ 


10-^ 


10-7 


10-7 


— > 9 


10-11 


10-11 


10-11 


io-« 


10-11 


10-11 


9 (quark) 




10-13 


10-11 


10-9 


10-13 


10-13 


9 (gluino) 


10-15 


10-15 


10-14 


10-11 


10-15 


10-15 
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